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Abstract 

We have studied the Giant Magnetoresistance (GMR) for laterally confined 

multilayers, e.g., layers of wires, using the classical Boltzmann equation in the 

current-in-plane (CIP) geometry. For spin-independent specularity factors at 

the sides of the wires we find that the GMR due to bulk and surface scattering 

decreases with lateral confinement. The length scale at which this occurs is 

of order the film thickness and the mean free paths. The precise prefactor 

depends on the relative importance of surface and bulk scattering anisotropics. 

For spin-dependent specularity factors at the sides of the wires the GMR can 

increase in some cases with decreasing width. The origin of the change in 

the GMR in both cases can be understood in terms of lateral confinement 

changing the effective mean free paths within the layers. 
PACS numbers: 75.70.-i, 75.70.Cn, 75.70.Pa, 72.15.Gd, 73.50.-h 
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I. INTRODUCTION 



Electrical transport properties of magnetic multilayers, which are thin alternating layers 
of ferromagnets (FM) and paramagnets (PM), has drawn considerable interest in recent 
years.0 A large decrease in the resistivity from antiparallel to parallel alignment of the film 
magnetizations has been observed experimentally.!'! This phenomenon is known as the giant 
magnetoresistance (GMR). The decrease in the resistivity from antiparallel to parallel align- 
ment arises from spin-dependent scattering.! The two sources of spin dependent scattering in 
these multilayers are bulk scattering and surface scattering. There have been numerous ex- 
perimental and theoretical studies to understand the physics of the giant magnetoresistance 
and to use it in applications.!'! 

There are several theoretical approaches to transport in magnetic multilayers. One ap- 
proach is to use the phenomenological Fuchs-Sondheimer theory of thin film resistance.!'! 
This approach is based on the classical Boltzmann equation. Other approaches are based 
on the linear response theory! - ! and the quantum Boltzmann equation.0 It has been shown 
that the conductivities obtained from the classical Fuchs-Sondheimer theory are in good 
agreement with the quantum results obtained via the Kubo formula.0 Thus it is useful to 
use the semiclassical approach to understand the physics of the giant magnetoresistance. 
Furthermore, this method with spin-dependent interface scattering^! reproduces the quali- 
tative features of the giant magnetoresistance seen in the experiments.^! 

With the development of nanotechnology, it is becoming important to understand the 
effect of lateral confinement on the GMR.Ill In this paper, we address the following questions: 
(i) Does the GMR increase or decrease with the reduction in width? (ii) If there is an 
increase or decrease in the GMR what causes it? (iii) What are the relevant length scales 
in the problem? We use the classical Boltzmann equation to answer the above questions. 

The layout of the rest of our paper is as follows. In Sec. II a detailed description of 
the model and our numerical procedure is given. In Sec. Ill we present the results of our 
calculation, and Section IV contains the conclusion. 
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II. THE MODEL 



The geometry of our problem is shown in Figure 1. We compute the current-in-plane 
conductivity (CIP) for a three- wire structure stacked along the y-direction. Both the current 
and the electric field are in the z direction. We have two identical ferromagnetic materials 
(FM) and a paramagnet (PM). For this three wire geometry we compute the conductivity in 
the case when the ferromagnet's magnetizations are parallel, of, and when they are antipar- 
allel, oaf- The magnetoresistance is defined as the ratio of the change in the conductivity 
from parallel to antiparallel alignments divided by the parallel conductivity: 

GMR=1-— . (1) 
a F 

To calculate the conductivities for different spin alignments we consider the classical 
steady state Boltzmann transport equation: 

vV r /(v,r)--E.V v /(v,r)=[^l^] , (2) 

\ / scatt 

where /(v, r) is the distribution function for electrons of mass m at position r with velocity 
v in presence of the electric field E. Because of the scattering term, Eq. (|2]) is complicated. 
For simplicity we consider the case where the scattering term is: 

(df(v,T)\ /(v,r)-</(v,r)> 



dt J t 

/ scatt 



(3) 



with (/(v, r)) being the spherical average of the distribution function, and r is the relaxation 
time. This is the simplest scattering term to represent elastic scattering. To solve Eq. @j 
with the right hand side of Eq. (|3|), we define g(v, r) to be the deviation of the distribution 
function from its equilibrium value: 

/(v,r) = /*(|v|) + 0(v,r). (4) 

We next make the ansatz that within linear response the spherical average of the distribution 
function is equal to the equilibrium distribution function, (/(v, r)) = / e9 (|v|), i.e., the 
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spherical average of g(v, r) is zero. This will be checked explicitly later. The linearized 
Boltzmann equation for a wire labeled by an index n is then given by 

v ■ Vrf„(v, r) - -E • V v /S(|v|) = (5) 

171 T ns 

where s denotes the spin of the electrons. When the electric field is zero, Eq. (|5]) becomes 
a homogeneous equation with the expected solution being g ns (y,r) = 0. For non-zero 
electric field, g ns (v,r) is proportional to (eE/m)V v /^(|v|) = eE • vdf^j(\v\)/de v . Thus, 
g ns {— v, r) = —g ns (v,r), irrespective of the boundary conditions, and the ansatz that the 
spherical average of <? ns (v, r) is zero is justified. 
The general solution of the above equation^ is 

^(v,r)=^ s (v,r B )e-l r - rB l^l v l 



+ (T( ,,.E-v(^)(l-e-l r - rB l/ T -l v l), 

(6) 



where r# is a point on the boundary or interface. Eq. (|6|) implies that to find the distribu- 
tion function at position r with velocity v, we proceed from r backwards along v until we 
reach a point r# at the boundary. The distribution function at the boundary, g ns (v, tb), 
is determined by the boundary conditions. We recover the usual bulk value if we go far 
away from the boundary points. Also we can see from Eq. (§) that the electrons lose their 
momentum as they diffuse into the medium and the characteristic length scale for this is 
just the mean free path, r ns vp, where Vp is the Fermi velocity. 

We now examine the boundary conditions. At each interface the electrons undergo either 
specular or diffuse scattering. For the n th interface, which is between wires n and n + 1, 
we define the probability of spin s electrons being diffusively scattered as (1 — S ns ), where 
S ns is the spin dependent specularity factor. The probabilities for being specularly reflected 
and transmitted are S ns R ns and S ns T ns , respectively. The sum of R ns and T ns is one. The 
angular dependence of the surface scattering parameter ,E1 the reflection coefficients, and the 
transmission coefficiently has been studied, but in our calculation we treat those as angle 



independent. With these definitions, the boundary conditions at the nth interface can be 
expressed as: 

0£*(v,a;,y n ) = S ns T ns g™_ ls (v,x,y n ) 

+ S ns R ns g™(-v, x, y n ), (7) 
g%?u(-v,x,y n ) = S ns T ns g™(-v,x,y n ) 

+ S ns R ns g™_ ls (v } x, y n ), (8) 

where y n is the position of the interface, and the superscripts out, in correspond to electrons 
going out from the boundary or coming in to the boundary. Similar equations are satisfied 
at the sides of the wires except that there is no transmission. 

We use an iterative procedure to compute the distribution functions, g^f{v, tb), which 
are non- uniform along the x and y directions of the interfaces and edges (see Fig. 1). 
From Eq. (|J), we observe that the fi^f (v, r)'s with different velocities projected along the z 
direction are decoupled. This allows us to discretize the g%f(v, r) according to cos 9 = v z /\v\ 
and solve each separately. For each cos#, g°"'(v,r) at the edges and the interfaces carry 
two more indices: an angle and a position r^. The g°"*(v,r)'s for different and are 
related via Eq. (|6]) and the boundary conditions given in Eqs. (0) and (H). To illustrate this 
we consider the relations between the distribution functions of the outgoing electrons at the 
first interface in Fig. 1: 

0, 0, Tj) = SuRug^icose, 2tt - 0, r<) 
x e -di/nsv F sme + S la T u g°f(co86, 0, r k )e- d2/T2sVFSlne , 

(9) 

where d\, d,2 are the path lengths projected onto the x—y plane from one boundary to another. 
Initial guesses for g°"'(cos 9, 0, r^) are taken from the nearby g°"*(cos 9 — <5(cos 9), 0, rj). The 
calculations are converged to within 1% with the total number of divisions for cos#, and 
chosen as A^ cose =100, iV^=500, and iVj=400, respectively. 
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Once the distribution functions at the boundaries are known, by using Eqs. and 
(§) the distribution function of the electrons with momentum v at any point r can be 
determined. We explicitly check that our ansatz is valid: (g ns (v,r)) = 0. The current 
density along the direction of the electric field for wire n and spin s is: 




The conductivity is obtained by averaging the current density over a cross-sectional area, 
A: 

1 3 f 

a = -^jHYl Jns{x,y)dxdy. (11) 

ajK n=U=U 

Finally, the giant magnetoresistance is obtained from Eq. (|]) using parallel and antiparallel 
alignments of the magnetizations in Eq. ( |TTD . 



III. NUMERICAL RESULTS 

We begin this section with the different input parameters of our problem. They are: 
the spin dependent mean free paths in the ferromagnet, L*^, the mean free path in the 
paramagnet, L PM , the spin-dependent transmission coefficients, T^, the specularity factors, 
S, and the thicknesses of the layers, t FM and t PM . We choose the mean free paths to be 
those determined in an experiment on a Co/Cu/Co structure!! L^=55A, L^= 10A, and 
L PM = 226A. The transmission coefficients are obtained from an average of the transmission 
coefficients calculated by Stiles:0 = 0.8 and = 0.4. We take the specularity factors 
at all interfaces and sides to be the same: S = 0.9. Finally, the thicknesses are chosen 
to be those of Cu when the Co layers are antiferromagnetically coupled at zero field.ll! 
For simplicity we take the ferromagnets to have the same thicknesses as the paramagnets: 

^FM = tpM- 

With the above parameters the giant magnetoresistance is due to both spin anisotropies 
in the bulk and surface scattering. It is useful to consider the limiting cases when the GMR is 



6 



due to only bulk scattering anisotropies or only surface scattering anisotropies. To do this, we 
consider two special cases: (a) when the transmission coefficients are equal: = = 0.6 
and (b) when the bulk mean free paths in the ferromagnet are equal: = = 17 A. In 
the following these are referred to as the (a) bulk scattering case and (b) surface scattering 
case. 

In Fig. 2 we plot the giant magnetoresistance as a function of width for the three cases: 
(a) GMR due to bulk scattering anisotropies, (b) GMR due to surface scattering anisotropies, 
and (c) GMR due to both bulk and surface scattering anisotropies. In all cases the giant 
magnetoresistance decreases as we reduce the width. To understand this decrease with 
lateral confinement we consider a slab, i.e., a wire with infinite width. Diffusive scattering 
at the sides of the wire reduces the conductivity and mean free paths of the wire compared 
to the slab. This reduction in the conductivity in going from a slab to a wire should be 
comparable to the reduction in going from a bulk system to a slab with a thickness equal to 
w. The conductivity of such a slab is given byi 

(ne 2 L\f 3L 
\ vtlvf j { 2w 

f 1 ,/ n \ n •? „ 1 — exp(— w/Ll cos^l) 1 

x / d(cos0 cos sin 2 0- — P \ '-± %-\, 12 

Jo [1 — S exp(— w/L\ cos0|J] J 

and the bulk conductivity is 

(ne 2 L\ 

oWk = , (13) 

\mvF J 

where L is the mean free path of the electrons. We obtain an effective mean free path, 



L e g, by replacing L in Eq. flT3|) with L e ^ such that atmik = o"siab- The effective mean free 
paths can be used in a multilayer calculation with infinite width wires. The results of such 
an effective mean free path multilayer calculation are plotted as the solid lines in Fig. 2. 
This approximation is in good agreement with the exact results (symbols), and hence we 
conclude that the reduction in the GMR is due to a decrease in the effective mean free path 
within the layers. 
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Although the GMR is reduced in all the cases shown in Fig. 2, the actual GMR vs. width 
curves are different. In particular, there are different length scales at which the GMR is re- 
duced. The width at which the GMR is half its infinite width value, GMR(wire) / GMR(slab) 
= 1/2, is defined as the half- width. In Fig. 3 we have rescaled the GMR vs. width curves 
shown in Fig. 2 by GMR(slab) and the half-width. All the points fall close to a single 
curve. This means that the half- width and the GMR for infinite width wires, GMR(slab), 
determine the GMR vs. width curves. We have also tried the same rescaling with other 
values of S, which is still the same for all sides and interfaces, and found the same curve. 

In Figs. 4(a) - 4(c) we have plotted the half-width as a function of the mean free paths 
in the ferromagnet and the thicknesses. The cases (a) - (c) refer to the same parameters as 
in Fig. 2. The ratio between the mean free paths in the ferromagnet is fixed to 
in Figs. 4(a) and (c), which is the same ratio used in Figs. 2(a) and (c). As for Fig. 4(b), 
the mean free paths are equal in the surface scattering case of Fig. 2(b): V 1 = LK For 
bulk scattering anisotropics (Fig. 4(a)) the half- width increases almost linearly with both 
the thickness and the mean free path. On the other hand, for surface scattering anisotropies 
(Fig. 4(b)) the half- width depends primarily on the thickness of the film and only weakly on 
the mean free paths in the ferromagnet. When both surface and bulk scattering anisotropies 
are present (Fig. 4(c)), the dependence of the half-width on the mean free path and thickness 
can be complicated. In the region where the half- width has a peak in Fig. 4(c), the GMR 
has a local minimum as a function of the mean free path. The origin of this local minimum 
comes from the near cancelation of the GMR from the bulk and surface contributions. The 
location of this minimum depends primarily on the transmission coefficients. 

Although we believe the generic behavior is that the GMR will decrease when the mul- 
tilayers are laterally confined, one can find parameters where the GMR actually increases 
with lateral confinement. The two ways we have found to do this are (i) have the sides of 
the wires introduce additional spin dependence in the scattering and (ii) have the sides of 
the wires selectively decrease the resistivity in the ferromagnet relative to the paramagnet. 
The idea behind both of these is again that laterally confining the multilayers reduces the 

8 



effective mean free paths within each layer. By changing the mean free paths by different 
amounts, one can tune the GMR. As an example in Fig. 5 we have plotted the GMR vs. 
width for two cases which differ only by the specularity factor at the sides of the FM. For the 
solid curve the specularity factors are 0.9 for both spin-up and spin-down electrons, while 
for the dashed curve the specularity factor at the sides for spin-up electrons is 0.9 and for 
spin-down electrons is 0.5. The GMR actually increases as one decreases the width of the 
sample because the mean free path for spin-down electrons decreases more rapidly than the 
mean free path for spin-up electrons. Such spin dependent specularity factors may occur 
naturally or be attainable by coating the sides of the multilayer. Allowing spin-dependent 
specularity factors at non-transmitting interfaces opens up the possibility that one can create 
a GMR device without ferromagnetic conductors, but only insulating ferromagnets which 
change the surface scattering. 



IV. CONCLUSION 

In this paper we have studied the effect of lateral confinement on the giant magnetoresis- 
tance. For spin independent specularity factors at the sides of the wires, the GMR decreases 
as one reduces the width of the wires. For this case we found that the GMR vs. width curves 
could be very nearly collapsed onto a single curve by rescaling the GMR by its infinite width 
value and the width by the half-width. The half-width depends on both the mean free paths 
and the thickness of the films. In the case of the GMR due solely to surface scattering, the 
thickness of the sample plays a dominant role in determining the half-width, while in the case 
of GMR due solely to bulk scattering, the half-width increases roughly linearly with both 
the mean free paths in the ferromagnetic wires and the thicknesses. The general case when 
both surface and bulk anisotropies are important can lead to more complex dependencies. 

The source of the decrease in the GMR is a reduction in the effective mean free path 
in the layers due to scattering off the sides of the wires. We showed this quantitatively by 
determining an effective mean free path within each wire and substituting it into a multilayer 
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calculation for films (infinite width wires). The results of the approximate solution and the 
exact solution agree quite well. 

For the case of spin-dependent mean free paths we find that there are parameter regimes 
where the GMR can increase as the sample is laterally confined. The origin of the increase is 
again changing of the effective mean free path within each layer as one decreases the width. 
With spin-dependent specularity factors one can change the ratio of the mean free paths 
for spin-up and spin-down electrons and hence change the GMR. Thus, with appropriately 
prepared sides of the wires, one may be able to increase or at least stem the decrease in 
the GMR. In any case the length scale at which the reduction in the GMR takes place is 
typically quite small, of order the mean free paths and the thicknesses of the layers. The 
effects discussed here will not be important until one goes to very small samples. 

This work was supported by DOD/AFOSR grant F49620-96-1-0026 and NSF frant 
DMR9357474, and the NHMFL. We thank J. Childress, T. S. Choy for helping us in different 
stages of our work. 
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FIGURES 




Width 

FIG. 1. Schematic diagram of the three wire structure. Wires 1 and 3 are the ferromagnets of 
thickness tpM, and wire 2 is a paramagnet of thickness tpM- Current is in plane (CIP) along the z 
direction, and the wires are stacked in the y direction. S represents the spin independent specularity 
factors at the sides of the ferromagnet, and the paramagnet. To determine the distribution function 
for outgoing electrons at point ij, g(cos8,(p,Tj), one must consider all possible incoming electrons 
from points on the edges and interfaces, e.g., points r, and r&. The precise relationship is described 
in the text (see Eq. (|9|)). 
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10 20 30 
Width (A) 

FIG. 2. Giant magnetoresistance of a three wire structure as a function of width for (a) bulk 
scattering, (b) surface scattering and (c) both bulk and surface scattering. The symbols refer to 
different film thicknesses, tpM = ^fm = 8A (circle), 20A (box), and 30A (triangle). In all cases the 
GMR decreases as we laterally confine the multilayers. The origin of this decrease in the GMR can 
be understood in terms of changing the effective mean free paths in the wires. As the wire width 
is reduced the effective mean free path within each wire decreases. To make this more quantitative 
we obtain an effective mean free path for each wire and use these mean free paths in a multilayer 
calculation (infinite width wire). The results, which are shown as the solid lines, are in good 
agreement with the exact calculation (symbols). For cases (a) and (c) the mean free paths are 
chosen for a Co/Cu/Co structure,^ which has L^=55A, L^=10A, and L PM =226A. For case (b), in 
which only surface scattering contributes to the GMR, we use LT = L l =nA and L PM =226A. For 
cases (b) and (c) the transmission coefficients at the interfaces are taken to be TT=0.8, T^=0.4 
while for case (a), where the GMR is due only to bulk scattering, we take 7^=7^=0.6. The 
thicknesses for wire 2 are chosen to be that of Cu when the Co slabs are antiferromagnetically 
coupledjl! and for simplicity we choose the Co layers to have the same thickness as the Cu. In all 
cases the sides and the interfaces have a specularity factor 5=0.9. 
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FIG. 3. Rescaled giant magnetoresistance as a function of width. For large width the GMR 
for a wire approaches that for an ordinary unconfined multilayer, GMR(slab). We define the 
width at which the GMR is reduced to half of GMR(slab) as the half-width. Rescaling the giant 
magnetoresistance by GMR(slab) and the width by the half-width, the GMR vs. width curves of 
Fig. 2 fall onto a single curve. The symbols are the same ones used in Figs. 2(a)-(c). 
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FIG. 4. Half-width as a function of the thickness and mean free path in the ferromagnet 
(L-I-) for (a) bulk scattering, (b) surface scattering, and (c) bulk and surface scattering. In Fig. 
4(a) and (c) the ferromagnetic spin-up and spin-down mean free paths are kept at a fixed ratio, 
L^/L^ = 5.5, which is the same ratio used in Figs. 2(a) and (c). In the surface scattering case 
of Fig. 4(b) ]J equals LK The half-width depends both on the thickness and the mean free path 
for all three cases. In case (a), bulk scattering, the half- width increases roughly linearly with the 
thickness and the mean free path, while in case (b), surface scattering, the half- width depends 
primarily on the thickness of the layers. The general case when both bulk and surface scattering 
are important (case (c)) can lead to complex dependence on the mean free path and film thickness. 
The feature at small in (c) is associated with the fact that the GMR has a local minimum in 
this region (see discussion in text). Except for t and L^, which vary, the parameters used in (a), 
(b), and (c) are the same as in Fig. 2. 
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FIG. 5. Giant magnetoresistance as a function of width for (a) spin-independent scattering 
at the sides and interfaces, (b) spin-dependent scattering at the sides of the ferromagnet. For 
(a), the GMR decreases with decreasing width (solid line), whereas for (b), the GMR increases 
with decreasing width (dashed line). The increase of the GMR with reduced width is due to the 
spin-down mean free path decreasing faster than the spin- up one. For (a), we choose the same 
parameters as of Fig. 2, except the mean free paths for the spin-up and down electrons in the 
ferromagnet are 200A and 100A, respectively. For (b), we use the same parameters except that 
the specularity factors for the sides of the ferromagnet are spin-dependent: 5j ide =0.9 and 5g idc =0.5. 
Note that the GMR eventually goes to zero for small enough widths because the effective mean 
free path in the paramagnet becomes much smaller than the thickness of the paramagnetic layer. 
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